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It is shown that current saturation in semiconducting carbon nanotubes is indistinguishable from 
metallic nanotubes if the carrier density is above a critical value determined by the bandgap and 
the LO phonon energy. This feature stems from the higher number of current- carrying states in the 
semiconducting tubes due to the van-Hove singularity at the band-edge. Above this critical carrier 
density, the ensemble saturation velocity at high- fields is found to be independent of the bandgap, 
but strongly dependent on the carrier density, explaining recent observations. 

PACS numbers: 81.10.Bk, 72.80.Ey 



Introduction: Carbon nanotubes (CNTs) exhibit ex- 
ceptional current carrying capability. For example, Yao 
et al. in a seminal work on single-wall metallic CNTs 
(m-CNT) observed that the current carrying capacity, 
as high as 10 9 A/cm 2 , was fundamentally limited by 
the emission of LO phonons by energetic electrons at 
high electric fields [1]. The exceptionally high current 
carrying capacity was attributed by the authors to the 
fact that each carrier within a energy bandwidth de- 
termined by the LO phonon energy Hlulo ~ 160 meV 
moves with the Fermi velocity (vp °° 10 8 cm/s), lead- 
ing to a fundamental limit of the current determined by 
I = (4e/L)Y, k fkVF = 4e/ LO , where f LO = ^lo/2tt, 
L is the CNT length, the factor g = g s g v = 4 is the 
product of the spin and valley degeneracy, and fa is the 
occupation probability of state \k). The current evalu- 
ates to « 25 /iA, precisely explaining the experimental 
saturation current. This result is fundamental for long 
m-CNTs, where transport is in the diffusive limit; for 
short m-CNTs, the current can exceed the diffusive limit 
of 25 /iA due to ultrafast phonon emission leading to non- 
uniform heating. It has been shown that the increase in 
the saturation current over this limit can be explained by 
consideration of non-equilibrium (hot) phonons [2, 3]. 

There have been reports of similar current saturation 
in semiconducting CNTs (s-CNT) as well [4, 5, 6]. How- 
ever, the physical mechanism of current saturation is not 
yet clear. Previous theoretical considerations indicate a 
velocity saturation [7], and in [8], a dependence of the 
ensemble saturation velocity on the carrier density was 
established by the solution of the Boltzmann Transport 
Equation (BTE). The purpose of this letter is fourfold: 
a) to show that for ID s-CNTs current saturation is as 
fundamental as in m-CNTs, b) as opposed to m-CNTs, 
the saturation current in a s-CNT is sensitive to the lo- 
cation of the Fermi level Ep before the application of a 
bias, c) the net ensemble velocity ('saturation velocity') 
of carriers in s-CNTs is not universal, but strongly depen- 
dent on £p as well, and can range from 0.1-lx^, and 
d) the saturation velocity of carriers in s-CNTs becomes 
independent of the bandgap for large carrier densities. 
These four facts are proved using analytical arguments 
based on the bandstructure, and the ultrafast electron- 
optical phonon interactions in CNTs. They offer a simple 



explanation for various recent observations made from 
experiments and numerical simulations [7]. 

CNT bandstructure & saturation current: For the anal- 
ysis of transport properties, it suffices to consider expan- 
sions of the energy dispersion for small momenta around 
the /C, KJ points of the underlying graphene bandstruc- 
ture [9] . Let the bandstructure of the n th subband of a s- 
CNT be given by the dispersion E n (k x ) = hvp \/^ + ^, 
where h is the reduced Planck's constant, k x is the longi- 
tudinal wavevector, and k n is the transverse wavevector 
quantized by the diameter of the CNT. All energies are 
measured with respect to the Dirac point of the under- 
lying graphene bandstructure. The bandgap of the CNT 
with this bandstructure is given by E g = 2hvpko, where 
fco is the smallest allowed (nonzero) transverse wavevec- 
tor, fixed by the size quantization (£ g ~ 800 jd meV, 
where d is the diameter of the s-CNT in nm). 

An electric field F = — Fqx is applied across the drain- 
source contacts to the s-CNT, which are assumed to be 
perfectly transparent to electron flow. To determine the 
exact shape of the new Fermi surface, one needs to solve 
the BTE, as has been done in [8] for T = 300 K. However, 
for the T — » limit, various analytical results may be 
derived to illustrate the physics of current saturation in 
s-CNTs. The Fermi surface at low temperatures is sharp, 
and the force due to the electric field (— e)F = hdk/dt 
populates +k x states by emptying out — k x states near 
the Fermi surface. If the contacts do not inject carri- 
ers into the wire (no space-charge currents), the equilib- 
rium occupation function = Q(kF — \k x \) (where 0(...) 
is the Heaviside unit-step function) shifts to the non- 
equilibrium position fa = 1 for — k^ < k x < while 
conserving the particle number inside the wire. Here 
kL,kR are the wavevectors of the highest right-going and 
the highest left-going electrons in response to the electric 
field. This is shown schematically in Fig 1. 

The density of states can be split equally into the left- 
going and right-going states: p to t(£) = Pl(£) + Pr(£), 
where 

Pr,l( £ ) = \ ' -4~, / = = > (!) 

2 TXhV F I g 2 _ (£jl)2 

and the resulting right- and left-going charge carrier 
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FIG. 1: Schematic representation of the mechanism of LO 
phonon scattering in s-CNTs. If the Fermi level before the 
application of a bias is high in the conduction band, the ultra- 
fast optical phonon scattering leads to the steady state carrier 
distribution shown in the middle and right figures; the satu- 
ration current is effectively carried by right-moving electrons 
spread over an energy bandwidth of Hulo- 




FIG. 2: The occupation function of carriers for a s-CNT with 
n < n cr at high fields. The emission of optical phonons is 
possible when the highest occupied state and the lowest unoc- 
cupied state differ by an energy Huolo] the lowest unoccupied 
state here is the band edge. The group velocity of carriers as 
a function of their energy is also shown. 



densities can be then written in the form ur^l 

gR,L 
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which defines the quasi- Fermi levels £j? and £p. Before 
the application of a bias, £p = £p = £ F , and the ID 
carrier density is given by 



no = n R + n L 
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If the contacts to the s-CNT do not inject excess carriers 
into the tube, then to ensure particle number conserva- 
tion no = nn + til requires that at all bias conditions, 
the relation 
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must be satisfied for the quasi- Fermi levels. 

At a high field, the distribution function reaches a 
stage when the difference between the highest filled elec- 
tronic state (HFES) and the lowest empty electronic state 
(LEES) equals the LO phonon energy to allow for energy 
relaxation by the emission of optical phonons. Building 
upon the theory in [1], we make a hypothesis that if the 
LO phonon emission process is ultrafast, then the dis- 
tribution function is locked in this configuration, and is 
resistant to any further increase in applied bias. Depend- 
ing on the availability of charge carriers, the LEES may 
be £p, or the band edge {£ = £ g /2). Before we proceed to 
calculate the net saturation current, we introduce some 
critical parameters that serve to highlight this fact. 

Critical parameters: A critical Fermi level £ F cr is the 
equilibrium Fermi energy such that at a high bias, £j? — 
£p = Tiwlo, and £j? = £ g /2. From Eq. 4 the critical 
Fermi energy must be 



f° - 



-F.cr ■= 2 \l (^LO) 2 + tkV LO £g + £% , 



(5) 



and is dependent only on the bandgap and the LO 
phonon energy of the s-CNT, which are fixed by the di- 
ameter and the lattices structure. The corresponding 
equilibrium ID critical carrier density is 



hcj LO (hu LO + £ g ), 



(6) 



which again is a fundamental quantity for the s-CNT. 
Physically, if the carrier density is above this, the LEES 
is a left-going state (as depicted in Figure 1), and if the 
carrier density is lower, then the LEES is is a right-going 
state, as depicted in Figure 2. 

We define two ratios: a = £ 9 /2£p and (3 = Hujlo /2£p. 
In addition, we consider only those semiconducting tubes 
for which £ g > Jilolo so that no interband optical 
phonon-assisted transitions are allowed. This condition 
implies that (3 < a < 1 when £p > £ g /2, i.e., when the 
Fermi level is inside the band and carriers are available 
for current conduction. If £p < £ g /2, no current flows at 
T^OK. 

Case 1: High carrier density: We first consider the 
situation when £p > £% cr , and n > n cr . Again, by 
ensuring particle number conservation and the condition 
£p—£p = huiLO, the high-bias quasi- Fermi levels for the 
right and left going states are found to be 
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Case 2: Low carrier density: If on the other hand 
the equilibrium Fermi level is such that £p < £p cr , and 
n < n cr , then the condition £ J? — £ p = Hulo cannot 
be satisfied due to insufficient carriers in the conduction 
band. However, since the LEES is now a right-going 
state, the condition £p — £ 9 /2 = Hujlo still allows for 
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FIG. 3: Saturation current of s-CNTs (8 g =0.3 & 0.4 eV 
as a function of 8%. The insets show (a) the critical Fermi 
level 8 FjCr and (b) the critical carrier concentration n cr as a 
function of the bandgap of the s-CNT. 
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FIG. 4: Saturation velocity of carriers in s-CNTs with 8 g = 
0.2, 0.3, & 0.4 eV as a function of the ID carrier density. 
For carrier densities no > n cr , v sa t — 4/lo/^o, and has no 
dependence on the bandgap of the s-CNT; for no < n cr , the 
saturation velocity has a weak dependence on 8 g . 



the emission of LO phonons. Under this situation, the 
quasi-Fermi levels are given by 



E$ = Eg/2 + hLj LO , 



(9) 



and 



= £% - 3a 2 + 4/3 2 + 4a/3 - 8vWl - a 2 y/a - /?, 

(10) 

In this condition, the HFES 8 F gets pinned at the 
energy 8 g /2 + Hujlo^ and the LEES is at the bottom of 
the conduction band. This case is depicted schematically 
in Figure 2. 

Saturation Current: The current flowing through the 
s-CNT is given in the T — > K limit by the relation 
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where L is the length of the CNT, and v g (k) = 
h~ x V k£(k) — vpk x j 'y/k% + k% is the projection of the 
group velocity of the carriers in the direction of the elec- 
tric field (Fig 2). The saturation current for the CNT 
evaluates to 

Isat(k n , k F ) = e±V F [yJk% + kl- yjkl + fc 2 ], (12) 

which may be reduced the simpler form 

r ^9^F 
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This is no different from what has been shown in 
[1] for m-CNTs, but it also holds for s-CNTs. The 
key point here is that for low fields, corresponding to 
zVds < h^LOi Ep — 8^ = cVdSi an d one recovers the 



Landauer relation / = g(e 2 /h)VDS- However, the same 
relation may be used to understand saturation currents 
at high bias conditions as well. Since the respective quasi- 
Fermi levels were shown earlier to depend on the equilib- 
rium Fermi level (or, indirectly, the carrier density), the 
saturation current in a s-CNT depends on 8 F . This is 
in stark contrast to m-CNTs, for which the quasi-Fermi 
level separation is Hcolo leading to a saturation current of 
I sat = tefho ~ 25/iA, no matter where the equilibrium 
Fermi level 8 F . 

In a s-CNT, the group velocity of a fraction of carriers 
near the band-edge is less than or equal to vp, whereas for 
a m-CNT is it always equal to vp. Thus, the single sub- 
band saturation current in a s-CNT can never exceed the 
current in a m-CNT for the same carrier density. For high 
ID carrier densities in s-CNTs (corresponding ton > n cr 
or equivalently 8 F > £ Fcr derived in Eqs 5, 6), the sat- 
uration current is identical to that in a m-CNT, i.e., 
I sa t = 4e/LO since the condition 8 F — 8 F = Tujlo holds. 
This similarity can be understood by noting that there is 
a van-Hove singularity at the band-edge of s-CNTs. This 
implies that there are many more carriers that contribute 
to the saturation current in a s-CNT than in a m-CNT, 
though a fraction of those carriers move with velocities 
less than the Fermi velocity. Since the DOS for a m-CNT 
is Pm(£) — g/^tiVF, the number of right-going carriers 
over a energy bandwidth Jilolo is um = ^fhol^F- The 
ratio of effective current-carrying states at saturation in 
a high-density s-CNT to that in a m-CNT is therefore 
greater than n cr /nM = yl + (£ g /^Lo)- 

For carrier densities in s-CNTs less than n cr , the satu- 
ration current is given by Eq. 13, with 8 F ,8 F from Eqs. 
9 & 10. This dependence exemplifies how a s-CNT un- 
der a high source-drain bias 'switches-ofT' when the Fermi 
level is pulled down towards the band-edge and then into 
the gap, for example, by electrostatic gating. The de- 
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FIG. 5: The velocity- field curves for carriers in s-CNTs of 
two bandgaps (solid line: £ g = 0.3 eV, dashed line: £ g = 0.4 
eV. The strong dependence of the velocity saturation on the 
effective carrier concentration, as measured by £% — 8c is 
indicated. The low- field mobility assumed is /io = 40, 000 
cm 2 /V.s for the calculation. 



pendence of the saturation current on £ F is plotted in 
Fig 3 for s-CNTs for two representative bandgaps. The 
insets show the critical Fermi level and the critical ID 
carrier concentration above which the saturation proper- 
ties of the s-CNT become indistinguishable from that of 
a m-CNT. 

Saturation velocity and velocity -field curves: Under the 
assumption that all carriers in the conducting band con- 
tribute to the saturation current equally, the saturation 
velocity of carriers may be written as v sat = I sa t/eno, 
which evaluates to 

— \(£f ~ £f) 9 (p r p L\ / 1zL \ 

Vsat - V F ■ == - • [t F - l F ), (14) 

^(£0)2 _ ( ^)2 ^hn 

and is strongly dependent on the ID carrier density in 
the band. Fig 4 shows the dependence of the saturation 
velocity on the density as well as the corresponding equi- 
librium Fermi level for s-CNTs of varying bandgaps. For 
carrier densities n < n cr , it increases and approaches v F 
as n — > and £ g decreases. However, from Eq. 14, the 
saturation velocity at carrier densities n > n cr is given 
by v sat = 4/lo/^o, and is inversely proportional to the 
ID carrier density. The theory is able to explain why the 
saturation velocity becomes independent of the bandgap 
of the s-CNT, as was found in numerical simulations of 



the BTE in [7]. 

If the low- field carrier mobility is /io, then a critical 
electric field may be defined as F cr = v sat /fi; the net 
velocity-field curve may then be approximated by the 
relation 

Vsat /-j r\ 

V ~ l + (v sat /^F)' [10> 

which is strongly dependent on the carrier density and 
the bandgap of the s-CNT. Fig 5 shows a few velocity- 
field curves plotted for s-CNTs of bandgaps 0.3 and 0.4 
eV, for different carrier densities, as indicated by the en- 
ergy £ F — £q. The low- field electron mobility chosen for 
this calculation is /io = 40,000 cm 2 /V.s, which is deter- 
mined primarily by elastic scattering processes (acoustic 
phonons & impurities). It is emphasized that the often- 
used concept of saturation velocity is an ensemble prop- 
erty, and care must be exercised in using it to predict 
carrier transit times, switching speeds, and allied quanti- 
ties since there are carriers that move at velocities higher 
as well as lower than this velocity. 

Deviations from the theory: Within the limits of the 
theory presented, the saturation current in a s-CNT due 
to single subband conduction in the diffusive limit can 
never exceed Io = ^c/lo- Any experimental observation 
of higher saturation currents in s-CNTs mustbe therefore 
attributed to effects that have not been considered. Five 
possible factors are listed: a) non-uniform heating due to 
ultrafast LO phonon emission [2] and/or hot-phonon ef- 
fects [10], b) ballistic effects (i.e., if the length of the CNT 
is shorter than the mean-free path for LO-phonon emis- 
sion L < Clo) [11] 7 c) occupation of multiple subbands, 
d) contribution from band-to-band Zener tunneling [12] 
and/or impact ionization, and e) space-charge injection 
into the s-CNT by the contacts. These factors can be in- 
corporated into the analytical framework presented here, 
but is not considered in this letter. 

Conclusions: In conclusion, a simple analytical theory 
is presented that shows that current saturation in s-CNTs 
and m-CNTs are indistinguishable if the carrier density 
in the s-CNT is above a critical value determined by the 
bandgap and the LO phonon energy of the s-CNT. The 
ensemble saturation velocity is found to be independent 
of the bandgap of the s-CNT for such carrier densities 
as well. The author acknowledges support from the NSF 
CAREER award (DMR-0645698), and useful discussions 
with Dr. Chagaan Bataar (ONR). 
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